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Lattice effects in a cerium based ciathrate compound CesPcboGes with a cubic Cr23C6-type struc- 
ture have been investigated by uitrasonic and thermai expansion measurements. Eiastic softenings 
of (Cii — Ci2)/2 and C44 proportional to the reciprocal temperature 1/T above Tqi = 1.25 K are 
well described in terms of the quadrupole susceptibility for the ground state Fg quartet. A huge 
softening of 50 % in (Cu — C\%)I1 and a spontaneous expansion AL/L — 1.9 x 10 -4 along the 
[001] direction in particular indicate the ferroquadrupole ordering of O2 below Tqi. The elastic 
anomalies associated with the antiferromagnetic ordering at Tn2 = 0.75 K and the incommensurate 
antiferromagnetic ordering are also found. Notable frequency dependence of C44 around 10 K is 
accounted for by the Debye-type dispersion indicating a Fs rattling motion of an off-center Ce ion 
along the [111] direction with eight fractionally occupied positions around the 4a site in a cage. 
The thermally activated T5 rattling motion obeying a relaxation time r = roexp(Fv/fc_gT) with an 
attempt time to = 3.1 x 10 -11 sec and an activation energy E = 70 K dies out with decreasing 
temperature, and then the off-center tunneling state of Ce ion in the 4a-site cage will appear at low 
temperatures. 

PACS numbers: 71.27.+a, 62.20.Dc, 65.40.De 



I. INTRODUCTION 

The 4f-electronic systems with spin and orbital de- 
grees of freedom in rare earth compounds frequently re- 
veal electric quadrupole orderings in addition to magnetic 
dipole orderings at low temperatures. The cubic com- 
pounds based on Ce 3+ ion with a Tg quartet ground state 
in particular have received much attention because the 
competitive phenomena associated with magnetic dipole, 
electric quadrupole and magnetic octupole degrees of 
freedom are expected. The direct product of r§ ® Tg 
is reduced to a direct sum Ti © T 2 © T 3 © 2T 4 © 2r 5 . 
The magnetic dipole J x , J y , J z belonging to T4 symme- 
try are order parameters for magnetic orderings. The 
quadrupole orderings of O®, 0\ with T 3 or O yz , O zx , 
O xy with T§ are interesting phenomena in the Tg sys- 
tem. We refer to CeAg exhibiting the ferroquadrupole 
(FQ) ordering of 0£ at T Q = 15 K. 1 - 2 CeB 6 is known 
as the antiferroquadrupole (AFQ) ordering of O xy -type 
with the propagation vector of fc=[lll] at Tq = 3.2 K. 3,4 
The octupole moments T xyz with T2 symmetry, T™, T™, 
T z with T 4 and Tg, , Tf with T 5 may play a role in 
the Tg system. 5 

A cerium-based ternary compound Ce3Pd2oGe6 with 
the Tg ground state has received much attention because 



the competition between quadrupole and magnetic or- 
derings is expected at low temperatures. 6 Ce3Pd2oGc 6 
crystallizes in a cubic C^Ce-type structure with a space 
group Fm3m consisting of four molecular units with 116 
atoms in a unit cell. 7 The twelve Ce sites located in 
cages are divided into two nonequivalent sites in crys- 
tallography. As shown in Fig. 1 the Ce ion at 4a site in 
a cage consisting of twelve Pd-atoms and six Ge atoms 
possesses point group symmetry Oh, while the Ce ion at 
8c site in a cage of sixteen Pd atoms has Td- The 4a 
sites form a face-centered cubic lattice, while the 8c sites 
make a simple cubic lattice. Inelastic neutron scattering 
on Ce3Pd2oGe6 revealed overlapping two peaks for the 
crystalline electric field (CEF) potentials, which corre- 
spond to magnetic dipole transitions from the Tg ground 
quartet to the Tr excited doublet at 60 K of the 4a site 
and from the Tg ground quartet to the T7 at 46 K of 
8c site. 8 The entropy obtained by low-temperature spe- 
cific heat measurement on Ce3Pd2oGe6 also indicates the 
ground state Tg quartet at both 4a and 8c sites. 9 

The low-temperature specific heat of Ce3Pd2oGe6 
shows a rounded small peak at Tqi = 1.25 K and a sharp 
A-peak at Tn2 = 0.75 K. 9 Magnetic susceptibility shows 
a clear cusp at Tn2, but exhibits no sign of anomaly at 
Tqi. 6 In addition to these experimental results, an elas- 
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FIG. 1: Cage of 4a site with Oh symmetry consisting of Ge 
and Pd2 atoms and cage of 8c site with T<j of Pdl and Pd2 
atoms in clathrate compound Ce3Pd2oGeg. The 4a site Ce2 
forms a simple cubic lattice, while 8c site Cel makes a face 
centered cubic one. The ferroquadrupole ordering below Tqi 
is relevant for the 8c sites and the antiferromagnetic ordering 
below Tn2 occurs at 4a sites. The Ts rattling motion origi- 
nates from the off-center Cel atom in 4a-site cage. 



tic softening of (Cn — Ci2)/2 in our preliminary paper 
suggests that the paramagnetic phase I transforms to the 
FQ phase II at Tqi and successively changes to the an- 
tiferromagnetic (AFM) phase III at Tn2- 10 The neutron 
scattering on Ce 3 Pd 2 oGe6 reveals a paramagnetic state 
of Ce ions at both 4a and 8c sites even in phase II be- 
tween Tqi and T N2 . The AFM ordering in phase III with 
a propagation vector k\ = [001] for cerium ions at 4a site 
is observed below Tn2- u Even in phase III below Tn2, the 
8c site still remains to be the paramagnetic state. The 
AFM ordering with incommensurate structure at 8c site 
appears only below T^ = 0.45 K. 

The clathrate compounds exhibiting the rattling mo- 
tion or off-center motion in a cage have received attention 
because their remarkable reduction of thermal conductiv- 
ity is favorable for application to thermoelectric device 
with a high figure of merit. 12 The ultrasonic waves are 
scattered by the rattling motion in an over-sized cage 
of a semiconductor Sr 8 Gai6Ge3o and a filled skutteru- 
dite compound PrOs4Sbi 2 . 13 ' 14 The off-center tunneling 
motion of OH ion doped in NaCl gives rise to elastic soft- 
ening at low temperatures. 15 The rattling motion in the 
present compound Ce 3 Pd2oGe 6 with clathrate structure 
has not been reported so far. 

In the present paper we show ultrasonic measurements 
on Ce3Pd 2 oGe6 in order to examine lattice effects asso- 
ciated with the quadrupole ordering and rattling motion 
in the system. The thermal expansion measurement is 
also employed to detect the spontaneous distortion be- 



low Tqi. In Sec. II, the experimental procedure and 
apparatus are described. The results of the elastic con- 
stant, magnetic phase diagram, thermal expansion arc 
presented in Sec. III. The ultrasonic dispersion due to 
rattling motion is also argued in Sec. III. In Sec. IV, we 
present concluding remarks. 

II. EXPERIMENT 

Single crystals of Ce3Pd2oGe6 used in the present mea- 
surements were grown by a Czochralski puling method. 
We have made the ultrasonic velocity measurements us- 
ing an apparatus consisting of a phase difference detec- 
tor. Piezoelectric plates of LiNb03 for the ultrasonic 
wave generation and detection are bonded on plane par- 
allel surfaces of sample. The x-cut plate of LiNb03 is 
available for transverse ultrasonic waves and the 36°j/-cut 
plate is for longitudinal waves. The ultrasonic velocity 
v was measured by fundamental frequencies of 10 MHz 
and overtone excitations of 30, 50 and 70 MHz. In the 
estimation of the elastic constant C = pv 2 , we use the 
mass density p = 10.254 g/cm 3 for Ce3Pd 2 oGc 6 with a 
lattice parameter a — 12.457 A. 6 

A homemade 3 He-refrigerator equipped with a super- 
conducting magnet was used for low-temperature mea- 
surements down to 450 mK in magnetic fields up to 12 
T. A 3 He- 4 He dilution refrigerator with a top-loading 
probe was used for the ultrasonic measurements in low- 
temperature region down to 20 mK in fields up to 16 T. 
Low input-power condition provides the low-temperature 
ultrasonic measurements free from a self-heating effect 
in the ultrasonic transducers. The sample length as a 
function of temperature or applied magnetic field was 
measured precisely by a capacitance dilatometer in the 
3 He- refrigerator. 

III. RESULTS AND DISCUSSIONS 

A. Temperature dependence 
of the elastic constants 

The elastic constants of Cn and Cl = (Cn + C12 + 
2C44)/2 of Ce3Pd 2 oGe6 in Fig. 2 were measured by the 
longitudinal ultrasonic waves with frequencies 10 MHz 
propagating along the [100] and [110] directions, respec- 
tively. The elastic constant (Cn — Ci2)/2 of Ce3Pd 2 oGc 6 
in Fig. 3 was measured by the transverse ultrasonic wave 
of 10 MHz propagating along the [110] direction polarized 
to the [110] one. The elastic constant C44 of Ce3Pd 2 oGe 6 
in Fig. 3 was determined by the transverse wave of 30 
MHz propagating along [100] polarized to [010]. The bulk 
modulus Cb = (Cn + 2Ci2)/3 in Fig. 2 was calculated 
by Cn in Fig. 2 and (Cn - C 12 )/2 in Fig. 3. 

It is remarkable that (Cn — Ci2)/2 exhibits a huge 
softening of 50% with decreasing temperature down to 
Tqi = 1.25 K. In phase II below Tqi the ultrasonic echo 
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signal of the (Cn — C\2)/2 mode completely disappears 
due to a marked ultrasonic attenuation. The softening 
of the longitudinal Cn and Cl modes in Fig. 2 origi- 
nates from the softening of (Cn — Ci2)/2, because Cn 
and Cl involve (Cn ~ C\2)/2 in part. The softening of 
(Cn — Ci2)/2 above Tqi and the spontaneous tetragonal 
distortion below Tqi, that will be shown in Sec. Ill D, 
provide evidence for the FQ ordering in phase II. The C44 
in Fig. 3 also exhibits a softening of 2.5% down to Tqi. 
The low-temperature behavior of Cn and C44 shown in 
insets of Figs. 2 and 3 indicates the transition to the FQ 
phase II at Tqi and successive transition to the AFM 
phase III at Tsr2 = 0.75 K. On the other hand, Cb shows 
monotonic increase with decreasing temperature. 

Neutron scattering on Ce3Pd2oGe6 revealed the para- 
magnetic state for Ce ions at both 4a and 8c sites in 
phase II, which is consistent with the present scenario of 
the FQ ordering at 8c site in phase II below Tqi. 11 The 
AFM ordering in phase III at 4a site below Tn2 has been 
detected by the neutron scattering. It has been proposed 
that the inter-site quadrupole interaction among 8c sites 
brings about the FQ ordering at 8c sites in phase II and 
Ce ions at 4a sites still remain to be the para-state even 
in phase II. The inter-site magnetic interaction among 
4a sites gives rise to the AFM ordering in phase III be- 
low Tn2- The magnetic ordering at 8c sites appears only 
below 0.4 K. We discuss about this transition in the fol- 
lowing Sec. Ill C. 



B. Quadrupole susceptibility 

In order to analyze the elastic softening of (Cn — 
C12V2 and C44 in Ce 3 Pd 2 oGe 6 of Fig. 3, we introduce 
the coupling of the quadrupole 0r 7 of Ce ions to the 
clastic strain £r 7 as 16 



Hq S = -^g r Or 7 (i)£r 7 , 



(1) 



where the summation ^T takes over Ce ions in unit 
volume and gr is a coupling constant. The inter-site 
quadrupole interaction mediated by phonons and con- 
duction electrons is written in a mean field approxima- 
tion as 
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FIG. 2: Temperature dependence of the elastic constants Cn, 
Cl and the bulk modulus Cb of Ce3Pd2oGe6. Inset shows 
the anomalies of Cn around the ferroquadrupole transition 
Tqi = 1.25 K and the antiferromagnetic transition Tn2 = 0.75 
K. 
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where (Or 7 ) denotes a mean field of the quadrupole 
and g' T means a coupling constant for the inter-site 
quadrupole interaction. By differentiating the total free 
energy consisting of 4f-electron and lattice parts with re- 
spect to the elastic strain er 7 , we obtain the temperature 
dependence of the elastic constant Cr (T) as 16 



Cr(T) = C° 



Nghr(T) 

i-g'rxATY 



(3) 



FIG. 3: Temperature dependence of the elastic constants 
(Cn — Ci2)/2 and C44 corresponding to shear waves in 
Ce3Pd2oGe6. Inset shows the anomalies of C44 around the 
ferroquadrupole transition Tqi = 1.25 K and the antiferro- 
magnetic transition Tn2 = 0.75 K. Solid lines are the calcula- 
tion by the quadrupole susceptibility for the T$ ground state 
and IV state at 46 K of Ce ions. The broken lines were back 
ground C° as shown in the text. A shoulder in C44 around 
30 K means the ultrasonic dispersion due to the T$ rattling 
motion. 
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Here Cp 1 (T) denotes a background elastic constant with- 
out the quadrupole-strain interaction and N is the num- 
ber of Ce ions in unit volume. The quadrupole suscepti- 
bility of xr(T) in Eq. (3) is written as 

-Ar(T) = (gg)-^{((gf>-(^W) 

where Ei is a second-order perturbation energy with re- 
spect to £r 7 for CEF state. 16 The first part in right hand 
of Eq. (4) corresponds to the Van Vleck-term and the 
second part to the Curie term. The Ce ions at both 4a 
and 8c sites in Ce3Pd2oGe6 have the T 8 ground state, 
while the state has excited energies of 46 K at 8c 
site and 60 K at the 4a site. As was already mentioned, 
the neutron scattering revealed that the FQ ordering of 
8c sites occurs at Tqi = 1.25 K and the AFM ordering 
of 4a sites appears at Tn2 = 0.75 K. These facts indi- 
cate that the inter-site quadrupole interaction of Eq. (2) 
among the Ce ions at 8c sites dominates the softening of 
(Cii— C\2)/2 as a precursor of the FQ ordering at Tqi. In 
the following analysis we simply assume the quadrupole- 
strain interaction of Eq. (1) and quadrupole interaction 
of Eq. (2) for the 8c site with the CEF splitting of T 8 (0 
K) and T 7 (46 K). 

The solid lines for (Cn — Ci 2 )/2 and C44 with Eq. 
(3) in Fig. 3 reproduce well the softening in param- 
agnetic phase I above Tqi. It should be noted that 
the softening above Tqi proportional to the recipro- 
cal temperature 1/T originates from the Curie term of 
Eq. (4). The coupling constants were determined to 
be |<7r 3 | = 107 K, g' r3 = 0.01 K for (Cn - C 12 )/2 
and |0r 5 | = 89 K, g' T5 = 0.19 K for C 44 . The back 
ground (C 1 ° 1 -C5 , 2 )/2 = (4.12 - 0.001T) x 10 10 J/m 3 and 
C° 4 = (3.44- 0.0007T) x 10 10 J/m 3 indicated by broken 
lines in Fig. 3 was used. The positive value of g r3 > 
are consistent with the FQ ordering in Ce3Pd2oGeg. A 
shoulder like anomaly in C44 around 10 K results from 
ultrasonic dispersion that is caused by a rattling motion 
of the rare-earth ion at 4a site in an oversized cage of 
Fig. 1. We discuss about this remarkable behavior in 
Sec. Ill E. 



C. Magnetic phase diagram 

In order to examine the magnetic phase diagrams of 
the FQ and AFM orderings in Ce3Pd2oGe6, we have 
made the low-temperature ultrasonic measurements of 
Cn, C44 and (Cn — C\2)/2 under magnetic fields. The 
softening of Cn of Fig. 4 reduces with increasing fields 
applied along the [001] direction parallel to the propa- 
gation direction of longitudinal wave. The FQ transi- 
tion points Tqi indicated by downward arrows in Fig. 4 
shift to higher temperatures and become indistinct in 
high fields up to 6 T. In Fig. 5, the FQ transition points 
Tqi also shift to higher temperatures accompanied by 
reduction of the softening in C44 with increasing applied 
fields along the [001] direction. 
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FIG. 4: Low-temperature behavior of Cn of CesPcboGee un- 
der magnetic fields along the [001] direction. Downward ar- 
rows indicate the ferroquadrupole ordering temperature Tqi 
and upward arrow means the antiferromagnetic ordering tem- 
perature Tn2- 
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FIG. 5: Low-temperature behavior of C44 of CesPd2oGe6 un- 
der magnetic fields along [001]. Downward arrows indicate the 
ferroquadrupole ordering temperature Tqi and upward arrow 
means the antiferromagnetic ordering temperature Tn2- 



In Figs. 6 and 7, we show the field dependence of 
(Cn - Ci 2 )/2 applying fields along the [001] and [110] 
directions, respectively. In zero magnetic field, the 
(Cn — Ci2)/2 mode exhibits the considerable softening of 
50% and the strong ultrasonic attenuation losing the echo 
signal in the vicinity of the FQ transition Tqi = 1.25 K. 
The magnitude of the softening decreases abruptly with 
increasing fields along both [001] and [110] directions. 
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FIG. 6: Low temperature behavior of (Cu — Ci2)/2 of 
Ce3PcboGe6 under magnetic fields along the [001] direction. 
The successive phase transitions I-II-III are indicated by ar- 
rows. 
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FIG. 7: Low temperature behavior of (Cu — Ci2)/2 of 
Ce3Pd2oGe6 under magnetic fields along the [110] direction. 
The successive phase transitions I-II-II' are indicated by ar- 
rows. 



In magnetic fields, clear minima of (Cu — Ci2)/2 corre- 
sponding to the transitions Tqi from the paramagnetic 
phase I to the FQ phase II shift to higher temperatures. 
Only broad round anomalies around Tqi have been ob- 
served in high fields of 5 T. This behavior of Tqi is similar 
to the results of the FQ ordering accompanied by the soft 
C44 mode in HoB 6 . 17 The Neel temperature Tn2 in Fig. 6 
shifts slightly to lower temperatures with increasing mag- 
netic fields. In Fig. 7 the anomalies associated with the 



transition between phase II and II' below Tqi have been 
found. 

For the investigation of low-temperature and high-field 
behavior in FQ II and AFM III phases, we have mea- 
sured field dependence of the C44 and C\ = (Cn + C12 + 
2Cj4)/2 employing the dilution refrigerator. In Fig. 8 we 
show C44 versus H at 30 mK in fields up to 12 T applied 
along [001]. Inset of Fig. 8 is expanded view below 2.5 T. 
An anomaly of the phase II-III boundary at 2.1 T indi- 
cated by a vertical line has been observed. Furthermore, 
several anomalies at 0.5, 1.2 and 1.6 T associated with 
sub-phases of the phase III have been found. It should be 
emphasized that appreciable hysteresis phenomena be- 
tween increasing and decreasing field sequences have been 
found only in phase III. 

In Fig. 9 we show the low-temperature field depen- 
dence of the Cl in fields up to 16 T applied along [110]. 
Low- field behavior below 2.5 T is shown in inset of Fig. 
9. We have observed a new phase boundary around 8.2 
T in phase II, which is probably a sub-phase of the FQ 
phase II. However, this phase boundary is absent in fields 
along [001] as shown in Fig. 8. At low field in phase III, 
we have found several anomalies in Cl of Fig. 9 showing 
a hysteresis behavior. As can been seen in inset of Fig. 9, 
this hysteresis becomes pronounced with decreasing tem- 
perature. These sub-phases with hysteresis behavior in 
phase III in magnetic fields along both [001] and [110] 
are well consistent with the results of neutron scattering 
experiments that detect weak incommensurate magnetic 
Bragg peaks with a propagation vector fe=[0 1-r], (r ~ 
0.06) at 8c site. 11 

The magnetic phase diagrams of Ce3Pd2oGe6 in 
Figs. 10 and 11 are obtained in fields along the [001] 
and [110] directions, respectively. We present the results 
of the ultrasonic measurements together with the results 
of thermal expansions in Sec. Ill D. It is of importance 
that the FQ phase II is stabilized in fields for the [001] 
direction in Fig. 10 and the [110] direction in Fig. 11. 
The FQ sub-phase II' was added to the phase diagram 
and the upper limit at 8.2 T of the phase II' newly de- 
termined in fields along [110] of Fig. 11. However, the 
FQ sub-phase is absent in fields along [001] of Fig. 10. 
This result indicates strong anisotropy of the quadrupolc 
interaction of 0% in Ce3Pd2oGe6- 

The series of R 3 Pd 2 oX 6 compounds usually show two 
successive AFM orderings of 8c site at Tni with a prop- 
agation vector k\ = [111] and of 4a site at Tn2(< Tni) 
with fc 2 = [001], Nd 3 Pd 2 oGe 6 {T m = 1.75 K, T N2 = 0.58 
K) 18 , Nd 3 Pd2 Si 6 (T N i = 2.4 K, T N2 = 0.7 K) 19 , 
Tb 3 Pd 20 Si 6 (T N i - 10.2 K, T N2 = 4.1 K) 20 , Dy 3 Pd 20 Si 6 
(Tni ~ 5.8 K, Tv 2 <~ 1.8 K) 21 and so on. One can rea- 
sonably expect that the transition temperature Tni at 
8c site is always higher than Tn2 at 4a site since the 
distance ~ 6.2 A between rare-earth ion of 8c site is 
much shorter than the one ~ 8.8 A of 4a site. In the 
present Ce 3 Pd 2 oGe6, at first the FQ ordering at 8c site 
with a structural change from cubic lattice to tetrago- 
nal one occurs at Tqi = 1.25 K. Therefore, the AFM 
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FIG. 8: Field dependence of the C44 at 30 mK in fields along 
[001] up to 12 T. Inset is expanded view below 2.5 T indicating 
the magnetic transitions. 
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FIG. 9: Field dependence of the C L = (C11 + C12 + 2C 44 )/2 
at various temperatures in fields along [110] up to 16 T. Inset 
is expanded view below 2.5 T indicating the magnetic transi- 
tions 



ordering at 8c site is hard to take place because the fa- 
vorable propagation vector k\ = [111] of 8c site docs 
not match to the tetragonal lattice in phase II. In other 
words, the AFM ordering at 8c site is replaced by the 
FQ ordering in Ce3Pd2oGeg. While, the AFM ordering 
of 4a site with a propagation vector k 2 = [001] is easy to 




T(K) 

FIG. 10: Magnetic phase diagram of CeaPd2oGe6 under fields 
along the [001] direction. The boundary from paramagnetic 
phase I to the ferroquadrupole phase II shifts to higher tem- 
peratures with increasing fields, while the boundary from 
phase II to the antiferromagnetic phase III shifts to lower 
temperatures in fields. 
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FIG. 11: Magnetic phase diagram of CeaPd2oGe6 under fields 
along the [110] direction. The sub-phase II' exists in the AFQ 
phase II. 



occur even in tetragonal structure below Tn2- Neutron 
experiments detected large enough value of saturation 
cerium moments /i(4o) = (1.1 ± 0.1)^B/Ce that is ex- 
pected from the ground state quartet T 8 perpendicular 
to the fei = [001] in Ce 3 Pd 20 Ge 6 far below T N2 at 50 
mK. 11 
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D. Thermal expansion 

In order to examine the structural change due to the 
FQ ordering at Tqi = 1.25 K, we have measured the ther- 
mal expansion along the [001] direction in Ce3Pd2oGe6- 
The sample lengths along the [001] and [111] directions 
are written by the symmetry strains as (Ai/i)r oi] = 
e zz = £b/3 + e u /V3 and (Ai/T)[ m] = e B /3 + 2(e yz + 
£zx + £xy)/3- Here, e B = e xx + e yy + e zz is a volume 
strain with Ti symmetry, e u = (2e zz — e xx — s yy )/V3 is 
a tetragonal strain with T 3 and e xy is a shear strain with 
T 5 . The length along [001] in Fig. 12 shows a monotonous 
decrease with decreasing temperature in paramagnetic 
phase I above Tqi and abruptly expands about AL/L = 
1.9 x 10~ 4 below Tqi. The thermal expansion along [001] 
in phase II and the huge softening of 50% in (Cn — Ci 2 )/2 
of Fig. 3 indicate the 2 ~tyP e FQ ordering accompanied 
by the structural transition from cubic lattice to tetrago- 
nal one with the spontaneous strain (e u ) in phase II. This 
spontaneous strain is proportional to the order param- 
eter as (e u ) = Ng r3 (0^)(2/(C^ 1 - C? 2 )) in mean-field 
approximation. Below Tn2 = 0.75 K, the AL/L along 
[001] slightly shrinks. Inset of Fig. 12 is expanded view 
of AL/L and the coefficient of the thermal expansion a 
at low temperatures. A sharp anomaly in the coefficient 
a has been found at the FQ transition Tqi. 

Measurements of AL/L versus T in various magnetic 
fields parallel to [001] are shown in Fig. 13. The magni- 
tude of the expansion AL /L in fields exhibits noticeable 
increase up to AL/L — 2.5 x 10~ 4 compared with that in 
zero magnetic field. The sharp increase of AL/L at the 



transition point to the FQ phase II has been observed in 
low magnetic fields below 1 T. On the other hand, the 
gradual increase in the thermal expansion AL/L above 
5 T up to 12 T indicates an obscure character of the I-II 
phase boundary in high fields. This is consistent with 
the fact that the clastic constants in fields of Figs. 4, 5, 6 
and 7 show obscure transitions in fields. This behavior is 
similar to the liquid-gas transition near the critical end 
point under hydrostatic pressures. 

The thermal expansion along the [001] direction in 
fields parallel to [001] and the considerable softening of 
(Cn - C 12 )/2 of 50% in Ce 3 Pd 2 oGe 6 strongly suggests 
that the order parameter of the FQ ordering in phase II is 
2 with T 3 symmetry. The relatively small softening of 
2.5% in C44 in Fig. 3 means that the quadrupole of O xy - 
type with T$ is irrelevant for the transition at Tqi. The 
thermal expansion of Ce3Pd 2 oGe6 along [111] is required 
to examine an interplay of the spontaneous strain e xy for 
the phase II. We refer to our recent study of the FQ tran- 
sition in H0B6 and the phase IV in Ce^Lai-xBe (x=0.75, 
0.70), 17 ' 22 where the trigonal strain e yz — e zx — e xy 
plays a significant role and the tetragonal strain e u is 
irrelevant. These facts are well consistent with the pro- 
nounced elastic softening in C44 of 70% in HoB 6 and of 
31% in Ce x Lai_ x B 6 (z=0.75, 0.70). 



E. Ultrasonic dispersion of C44 

The C44 mode associated with the elastic strain 
£yz , £zx , £xy of T 5 symmetry of Ce3Pd 2 oGe6 in Figs. 3 and 
5 exhibits a shoulder like anomaly around 10 K in addi- 
tion to the characteristic softening due to the quadrupo- 
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pansion AL/L at low temperatures. 



FIG. 13: Thermal expansion AL/L along [001] in 
Ce3Pd2oGe6 under fields up to 12 T along the [001] direc- 
tion . 
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lar coupling above Tq± — 1.25 K. It should be noted 
that this anomaly is absent for the (Cn — Cy^jl mode 
associated with r 3 elastic strain s v and the bulk modu- 
lus Cb with Ti volume strain £b- In order to examine 
the origin of this anomaly, we have measured the fre- 
quency dependence of 644 from 10 MHz up to 250 MHz. 
The elastic constant 644 of Fig. 14(a) exhibits shoulders 
showing remarkable frequency dependence. An increase 
in ultrasonic attenuation around shoulder has also been 
found, but not discussed here. We describe this frequency 
dependence of the elastic constant C^ 4 (lu) in terms of 
Debye-type dispersion as 

r<D 1 ,\ riDt nn \ ^44(00) - 6^(0) 

644(0;) = 6 44 (00) 1,22 ' 5 

1 + LO A T A 

where 6*44(00) and 644(0) are the elastic constants of 
high frequency limit and low frequency one, respectively. 
Here to is an angular frequency of the ultrasonic wave and 
t means the relaxation time of the system. In fittings of 
Fig. 14(b), we take into account the superposition of two 
susceptibilities by the quadrupole one of Eq.(3) and the 
Debye-type dispersion of Eq.(5) as 644 = Cf 4 + C^(w). 
The inflection points around 10 K in 644 indicated by ar- 
rows in Fig. 14(a) mean the temperatures where the r co- 
incides with the lo as lot — 1. The ultrasonic attenuation 
is expected to be maximum at the temperatures for lot = 
1. The solid lines of Fig. 14(b) being the calculations 
with Eq. (5) well reproduce the experimental results of 
Fig. 14(a). The relaxation time obeying the Arrhenius- 
type temperature dependence r = roexp(£ , /fcsT) with 
the attempt time tq = 3.1 x 10~ n sec and the activation 
energy E = 70 K has been determined. The parameter 
of AC = Cgioo) - 64^(0) = 0.004 x 10 10 J/m 3 is used. 

The ultrasonic dispersion due to electron thermal 
hopping has already found in the inhomogcneous va- 
lence fluctuation compounds of SUL3X4 (X=Se, Te), 
Yb4(As .7iSbo.29)3 23 and Sri 2 Ca2Cu2404i 24 . It is re- 
markable that the very slow relaxation time r and ex- 
tremely low activation energy E for Ce3Pd2oGe6 are ex- 
ceptional as compared to those of charge fluctuation com- 
pounds Sm 3 Se 4 25 and Sm 3 Te 4 26 , t ~ 2.5 x 10~ 13 sec and 
E ~ 1600 K, and Sr 12 Ca 2 Cu240 4 i, r = 1.01 x 10~ 13 sec 
and E = 1900 K. 24 This discrepancy of the order of tq 
and E between the present Ce 3 Pd 2 oGe 6 and the charge 
fluctuation compounds indicates that thermally activated 
rattling motion of heavy mass particle, which is proba- 
bly rare-earth ion in a cage, gives rise to the ultrasonic 
dispersion in Ce 3 Pd 2 oGe6. 

Glass materials and charge fluctuation compounds ex- 
hibit frequently the ultrasonic dispersion, which results 
from a thermally activated motion in a double- or multi- 
well potential. The compounds such as Sm 3 X 4 (X=S, 
Se, Tc) with different valence of Sm 2+ and Sm 3+ ions 
and Sri2Ca2Cu2404i with Cu 2+ and Cu 3+ ions cause the 
ultrasonic dispersions due to thermally assisted charge 
fluctuation in the temperature region between 100-200 K. 
The two-level system (TLS) due to an atomic tunneling 27 
or electron tunneling 26 manifests itself in glass materials 
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FIG. 14: Frequency dependence of the elastic constant of the 
transverse C44 mode in CeaPd2oGe6. (a) represents the ex- 
perimental results with the frequencies from 10 MHz up to 
250 MHz. (b) is the calculation in terms of the Debye-type 
dispersion of Eq. (5) in text. 



at low temperatures, where the thermally activated mo- 
tion dies out. The TLS yields the decrease in the elastic 
constant proportional to InT, 28 the specific heat to T 
and thermal conductivity to T 2 . 29 Besides in the case 
of Sm 3 Te4, remarkable logarithmic decrease in the elas- 
tic constant appears below about 15 K down to the spin 
glass transition at T g = 1.5 K, which suggests the exis- 
tence of the 4f-electron tunneling motion between Sm 2+ 
and Sm 3+ ions situated charge glass state. 26 

The present clathrate compound Ce 3 Pd2oGe6 is a crys- 
tal possessing an ideal periodic arrangement of cages 
in space. The stable trivalent Ce ions in cages of 
Ce 3 Pd2oGe6 are free from the valence fluctuation phe- 
nomena. As shown in Fig. 1, the clathrate compound 
Ce 3 Pd2oGe 6 is made up of the cage at 4a site con- 
sisting of Pd and Ge with distances dc c i-Gc=3.332 A, 
dcci-Pd2=3.067 A, and the cage at 8c site of Pd with 
dcc2-Pdi=2.868 A, dce2-Pd2=3.373 A. The trivalent Ce- 
ion with radii a = 1.7 ~ 1.8 A inside the 4a-site cage in 
particular is expected to show the rattling motion over 
the off-center positions being away from the center of the 
cage. Actually the neutron scattering on Pr 3 Pd 2 oGe6 
and Nd 3 Pd2oGe6 revealed the sharp transition peaks in- 
dicating the stable CEF splitting at 8c site and no indi- 
cation for CEF state at 4a site. 8 These results may imply 
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TABLE I: Rattling modes with eight off-center positions along the three-fold [111] direction in the 4a-site cage with Oh point 
group symmetry. Corresponding elastic strains are also listed. There is a Ts rattling mode coupled to the e^-type strain of 
the C44 mode. This is contrary to the absence of the T3 rattling mode to the e u ,e v of the (Cu — Ci2)/2 mode. The present 
ultrasonic dispersion in C44 of Cc3Pd2oGe6 originates from the IVtype rattling at 4a site. 



Symmetry 



Rattling mode 



Strain 



ri 
r 2 
r 3 

r 4 



pn — pi + P'2 + p3 + Pi + P5 + pa + P7 + ps 
PV'l — Pi + P'2 + p3 + Pi — PS — Pa — Pi — PS 



= Pi — Pi — Ps + Pi + Ph — Pa + P7 — Pa 
„ = Pi — P2 + ps - Pi + Ph - pa - pi + ps 
pr 4 ,z = pi + P2 — P3 — Pi — Pb — Pa + P7 + Ps 
Pr 5 ,yz = pi — P2 - P3 + Pi — Ph + Pa - pi + Ps 

n„ — n, — n., J- n„ — n. — n- -X- r,,. -I- n~ — n„ 



pn,i 
pn, y = pi 



Pr 5 ,yz = pi 

= pi — P2 + P3 — Pi — P5 + Pa + Pi — Ps 
pi + P2 — P3 — Pi + Ps + Pa — Pi — Ps 



PT 5 ,za 

Pr 5 



E~B — Exx Eyy Ezz 

£u = (2s 22 Exx ^t/t/)/V^3 
£v — — Exx ^yy 



E zx 

Exy 



the obscure CEF state due to the off-center Ce ion at 4a 
site contrary to the well-defined CEF splitting at 8c site 
being stable Ce ion position. 



F. T5 rattling motion 

Notable finding of Fig. 14 is that the ultrasonic dis- 
persion in the C44 mode associated with s xy -type strain 
indicates the rattling motion with specific T$ symmetry 
in Ce3Pd2oGe6. It is expected that the Ce ion in 4a cage 
with cubic symmetry Oh favors off-center positions along 
one of the three principle directions of [100], [110] and 
[111]. As one can see the cage at 4a site in Fig. 1, it is of 
particular interest that no atom exists along the three- 
fold [111] directions, while the Ge atom occupies along 
the four-fold [100] and the Pd atom along the two-fold 
[110] ones. This crystallographic character may promise 
a flat potential along the three- fold [111] directions and 
profound potentials along the four-fold [100] and two-fold 
[110] directions. Presumably the Ce ion at 4a site prefers 
the off-center eight positions along the [111] directions, 
which are defined as ri = (a, a, a), r 2 = (—a,— a, a), 
r 3 = {-a, a, -a), r 4 = (a, -a, -a), r 5 = (a, a, -a), 
r 6 = (-a, -a, -a), r 7 = (a, -a, a), r 8 = (-a, a, a). 
The atomic densities pi = pi(ri)(i = 1 -~ 8) at the 
eight off-center positions are also defined. When 48 sym- 
metry operators of Oh point group are acted on the 
atomic densities pi, one can derive the transfer repre- 
sentation matrices with 8x8 elements. Consequently, 
one obtains the characters X™/' f° r the rattling mo- 
tion by tracing the diagonal elements of the representa- 
tion matrices. Using the characters X™/' an d the char- 
acter table for the irreducible representations of Oh, 30 
the rattling motion over the eight off-center positions 
along the [111] direction is reduced to the direct sum 

of r 1 (iD)er 2 (iD)er 4 (3D)er 5 (3D). 

Employing projection operators on appropriate atomic 
density pi, we obtain the rattling modes for the ir- 
reducible representations as listed in Table 1 together 



with the elastic strains £r 7 - One can see the presence 
of T 5 rattling mode pn, V z, Pvb.zx, Pvb,x V coupled to the 
strain e yzi s zx , e xy contrary to the absence of T 3 rat- 
tling mode coupled to the strain e Ul e v . This is consis- 
tent with the fact that the ultrasonic dispersion reveals 
in C44 and is absent in (Cu — C\2)j2. The Ce atom 
in the cage of Ce3Pd2oGe6 obeys a harmonic oscillation 
of ((z) — (l/7rz )l/2exp(— z 2 /2zq) with a mean square 
displacement zq = (1/27t)(/ito/M) 1 / 2 . The attempt time 
t = 3.1 x 10~ n sec and the mass M — 140m p where m p 
is a proton mass, leads to the mean square displacement 
zq being approximately twice of off-center distances a as 
z Q = 2a = 0.48 A. 31 

The full symmetry Ti rattling mode pn = pi + pi + 
P3 + Pi + pb + Pe + p7 + Ps means the uniform atomic 
distribution with fraction 1/8 at each eight off-center 
positions. While the T 5 rattling mode, for instance 
PTb.xy = pi + P2 - Pa - Pi + P5 + Pe - Pi - Pd,, repre- 
sents anisotropic atomic distribution being quadrupolc 
O xy at the lowest order such as fraction 1 /4 at r\ , r 2 , r 5 , 
r 6 and zero at r 3 , r 7 , r 8 as shown in Fig. 15. The 
present group theoretical analysis for the rattling mode 
is essentially the same treatment previously argued for 
the charge fluctuation mode. 24 ' 26 ' 32 In the present system 
Ce3Pd2oGe6, the thermally activated T 5 rattling mode 
may be a ground state and the Ti, T 2 , T 4 be excited 
states. 

Our group has recently found similar ultrasonic dis- 
persion around 30 K in a heavy fermion superconductor 
PrOs4Sbi2 with a filled skuttcrudite structure. It should 
be noted that ultrasonic dispersion in the (Cu — Ci2)/2 
mode of PrOs4Sbi2 indicates the T 3 rattling motion of 
Pr atom over six fractionally occupied positions along 
[100]. The dispersion of (Cu - Ci 2 )/2 in PrOs 4 Sbi 2 is 
contrary to the one of C44 in the present compound of 
Ce 3 Pd 20 Ge 6 . 14 The attempted time t = 8.8 x 10~ n sec, 
activation energy E = 168 K and mean square displace- 
ment zo = 0.80 A in PrOs4Sbi2 are comparable to the 
present results of Ce3Pd2oGe6- 

The thermally activated T$ rattling motion with frac- 
tional atomic state in Ce3Pd 2 oGe6 dies out with decreas- 
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Pr5,xy=Pl+P2-P3-P4+Pr , -p6-P7-P8 

FIG. 15: Schematic view for the Fs rattling mode prs,x y due 
to the off-center Ce-ion along the three- fold [111] direction in 
the 4a-site cage. The prs,xy represents that fractional atomic 
density 1/4 is located at n, ri, re,, r$ and null at vs, r±, rj, 
rs- The freezing of the thermally activated motion of the Fs 
rattling mode brings about the atomic tunneling state at low 
temperatures. 

ing temperature. At further low temperatures, the off- 
center tunneling state of Ce ions in the 4a-site cages will 
appear, which means a quantum state being occupied 
four positions, for instance at n, r2, r$, for the case 
of pvb,xy, at the same time. The charge glass state in 
the inhomogeneous mixed valence compound Sni3Te4 re- 
vealed the low-temperature softening in elastic constants 
proportional to InT being resemble the structural glass. 26 
The present scenario of the atomic tunneling of Ce ions 
in 4a-site cages in Ce3Pd 2 oGe6 is also expected to show 
the elastic softening at low temperatures. However, the 
low-temperature quadrupole and magnetic orderings in 
Ce3Pd2oGe6 may mask the character of the tunneling 
state in the present case. In order to clarify the tun- 
neling and rattling in the present clathrate compounds, 
the low-temperature thermodynamic and ultrasonic mea- 
surements on La3Pd2oGe6 without 4f-electron in partic- 
ular is required. 

IV. CONCLUDING REMARKS 

In the present paper we have measured the elastic con- 
stants and thermal expansion of Ce3Pd2oGe6- The char- 
acteristic elastic softening in (Cn — Ci 2 )/2 and C44 is 
well described in terms of the quadrupole susceptibility 
for the r 8 ground state. The important finding is that 
the (Cn — Ci 2 )/2 shows the huge softening of 50% to- 
wards Tqi = 1.25 K and the C44 exhibits the softening of 
2.5% only. This result strongly indicates the FQ ordering 



with the order parameter of the T3 symmetry in phase II 
below Tqi. Actually we have successfully observed the 
sharp increase of AL/L = 1.9 x 10~ 4 in length along 
the [001] direction below Tqi. This is the evidence for 
the OSj'type FQ ordering accompanied by the structural 
change from cubic lattice to tetragonal one at Tqi in 
Ce 3 Pd 2 oGe 6 . 

For the investigation of the magnetic phase diagram 
concerning the FQ phase II and AFM phase III in 
Ce3Pd2oGe6, the elastic constants and thermal expan- 
sion in fields have been measured. We have found that 
the boundary from the paramagnetic phase I to the FQ 
phase II of O® shifts to higher temperatures with increas- 
ing magnetic fields. The result that the I-II phase tran- 
sition becomes obscure in fields is similar to the liquid- 
gas transition approaching to the critical end point under 
pressure. This result consistent with the fact that the FQ 
order parameter 0° in Ce3Pd2oGe6 has the total symme- 
try under fields along the [001] direction. 5 The boundary 
from the phase II to the AFM ordering shifts to lower 
temperatures as similar as the conventional AFM order- 
ing. 

We have found the ultrasonic dispersion in the C44 
mode indicating the rattling motion with T 5 symmetry. 
Taking into account the absence of the atom along the 
[111] direction in 4a-site cage, we have successfully picked 
up the specific T 5 rattling mode prb,yz, Pn.zx, Pn.xy with 
the fractional atomic density over the eight minimum 
positions of potential along the four- fold [111], [111], 
[111] and [111] directions. The dispersion of the C44 
mode obeying the Debye-formula revealed the thermally 
activated-type relaxation time r = Toexp(EfkBT) for the 
r$ rattling mode with an attempt time To = 3.1 x 10~ n 
sec and an activation energy E = 70 K. The estimated 
mean square displacement z — 0.48 A for the harmonic 
oscillation of Ce atom leads to the distance of the poten- 
tial minima along the [111] direction as a = zq/2 = 0.24 
A. In order to confirm the anisotropic atomic distribution 
in the 4a-site cage, the neutron or x-ray scattering is re- 
quired. The freezing of the thermally activated motion 
due to the T$ rattling mode with lowering temperature 
brings about the atomic tunneling state. By analogy of 
the charge glass compound Sm 3 Te4, the Ce-ion tunneling 
is expected at low temperatures. The ultrasonic investi- 
gation on La3Pd2oGe6 free from the long-range ordering 
due to 4f-electrons is now in progress to shed light on the 
tunneling and rattling in cages. 
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